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Abstract 

The emergence of the cosmological composition (the reheating era) after the inflationary period 

is analyzed in the framework of the braneworld models, in which our Universe is a three-brane 
00 ' 

. embedded in a five-dimensional bulk, by assuming the possibility of the brane-bulk energy exchange. 



o 



The inflaton field is assumed to decay into normal matter only, while the dark matter is injected into 
the brane from the bulk. To describe the reheating process we adopt a phenomenological approach. 



, by describing the decay of the infiaton field by a friction term proportional to the energy density of 

^Sj , the field. After the radiation dominated epoch the model reduces to the standard four dimensional 
> 

CS| ' cosmological model. The modified field equations are analyzed analytically and numerically in 

' both the extra-dimensions dominate reheating phase (when the quadratic terms in energy density 

l/-^ . dominate the dynamics), and in the general case. The evolution profiles of the matter, of the scalar 

o . 

OO , field and of the scale factor of the universe are obtained for different values of the parameters of 



the model, and of the equations of state of the normal and dark matter, respectively. The equation 



^ ■ describing the time evolution of the ratio of the energy density of the dark and of the normal 

H : 

matter is also obtained. The ratio depends on the rate of the energy fiow between the bulk and 
the brane. The observational constraint of an approximately constant ratio of the dark and of the 
baryonic matter requires that the dark matter must be non-relativistic (cold). The model predicts 
a reheating temperature of the order of 3 x 10^ GeV, a brane tension of the order of 10^^ GeV^, 
and the obtained composition of the universe is consistent with the observational data. 
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I. INTRODUCTION 



Early studies on superstring theory and M-theory have suggested that our four dimen- 
sional world is embedded into a higher dimensional spacetime. Particularly, the 10 dimen- 
sional (S) Eg heterotic superstring theory is a low-energy limit of the 11 dimensional 
supergravity, under the compactification scheme M^^ x S1/Z2 . Thus, the 10 dimen- 
sional spacetime is compactified as x CY^ x S1/Z2, implying that our universe (a brane) 
is embedded into a higher dimensional bulk. In this paradigm, the standard model particles 
are open strings, confined on the braneworld, whilst the gravitons and the closed strings can 
freely propagate into the bulk [3|. 

Among the collection of the braneworld scenarios, the Randall- Sundrum Type II model 
has the virtue of providing a new type of compactification of gravity . Standard 4D 
gravity can be recovered in the low-energy limit of the model, with a 3-brane of positive 
tension embedded in 5D anti-de Sitter bulk. The covariant formulation of the brane world 
models has been formulated in leading to the modification of the standard Friedmann 
equations on the brane. It turns out that the dynamics of the early universe is altered by 
the quadratic terms in the energy density and by the contribution of the components of the 
bulk Weyl tensor, which both give a contribution in the energy momentum tensor. This 
implies a modification of the basic equations describing the cosmological and astrophysical 

The recent observations of the CMB anisotropy by WMAP [8| have provided convincing 
evidence for the inflationary paradigm [9i], according to which in its very early stages the 
universe experienced an accelerated (de Sitter) expansionary phase (for recent reviews on 
inflation see 

At the end of inflation, the universe is in a cold and low-entropy phase, which is utterly 
different from the present hot high-entropy universe. Therefore the universe should be 
reheated, or defrosted, to a high enough temperature, in order to recover the standard Hot 
Big Bang HI]. The reheating process may be envisioned as follows: the energy density in 
zero-momentum mode of the scalar field decays into normal particles with decay rate F. The 
decay products then scattered and thermalize to form a plasma 10 1. 

Apart from the behavior of the inflaton field, the evolutions of dark energy and dark 
matter in reheating stage were also considered. In 121], dark energy and dark matter were 
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originated from a scalar field in different stages of the inflation, according to a special form 
of potential. Meanwhile, the conditions for unifying the description of inflation, dark matter 
a„d dark energy were eous.dered In Q. A spee.fie mode, was later proposed in Q. bv 
using a modified quadratic scalar potential. The candidates of dark matter in pj| and [141 ] 
were oscillations of a scalar field. However, it may be possible that dark matter existed on 
its own without originating from the scalar field. This may pose less stringent constraint on 
the scalar field, so that dark matter can be included in inflation paradigm in a easier way. 
On the other hand, it was proposed that the decay products of scalar field acquired thermal 
mass 

The reheating in the braneworld models has also been considered recently. In the context 
of the braneworld inflation driven by a bulk scalar field, the energy dissipation from the 
bulk scalar field into the matter on the brane was studied in 16|. The obtained results 



supports the idea that the brane inflation model, caused by a bulk scalar field, may be a 
viable alternative scenario of the early universe. The inflation and reheating in a brane 
world model derived from Type IIA string theory was studied in p/7|. In this model the 
inflaton can decay into scalar and spinor particles, thus reheating the universe. A model in 
which high energy brane corrections allow a single scalar field to describe inflation at early 
epochs and quintessence at late times was discussed in [l^. The reheating mechanism in 
the model originates from Born-Infeld matter, whose energy density mimics cosmological 
constant at very early times and manifests itself as radiation subsequently. The particle 
production at the collision of two domain walls in a 5-dimensional Minkowski spacetime was 
studied in 19(]. This may provide the reheating mechanism of an ekpyrotic (or cyclic) brane 
universe, in which two BPS branes collide and evolve into a hot big bang universe. The 
reheating temperature T^h in models in which the universe exits reheating at temperatures 
in the MeV regime was studied in |20|, and a minimum bound on Tjih was obtained. The 
derived lower bound on the reheating temperature also leads to very stringent bounds on the 
compactification scale in models with n large extra dimensions. The dark matter problem 
in the Randall- Sundrum type II brane world scenario was discussed in 21[, by assuming 
that the lightest supersymmetric particle is the axino. The axinos can play the role of cold 
dark matter in the universe, due to the higher reheating temperatures in the brane world 
model, as compared to the conventional four-dimensional cosmology. The impact of the 
non-conventional brane cosmology on the relic abundance of non-relativistic stable particles 
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in high and low reheating scenarios was investigated in 22|. In the case of high reheating 
temperatures, the brane cosmology may enhance the dark matter relic density by many order 
of magnitudes, and a stringent lower bound on the five dimensional scale may be obtained. 
In the non-equilibrium case, the resulting relic density is very small. The curvaton dynamics 
in brane- world cosmologies was studied in [23]. By assuming that the inflaton field survives 
without decay after the end of inflation, the curvaton reheating mechanism was applied to 
the Randall- Sundrum model and to its curvature corrections. 

Recently, motivated by particle physics results, the possibility of an energy exchange 
between the bulk and the brane was also suggested. The observational constraints on the 
cosmological models in the brane-world scenario in which the bulk is not empty, and ex- 



2|. The 



change of mass-energy between the bulk and the bane is allowed, were explored in 
evolution of matter fields on the brane is modified due to new terms in the energy momen- 
tum tensor describing this exchange. This model accounts for the observed suppression of 
the CMB power spectrum at low multipoles, and the observed recent cosmic acceleration is 
attributable to the flow of matter from the bulk to the brane. The cosmological evolution of 
a brane with a general bulk matter content was considered in {25! , by assuming that the bulk 
pressure and the energy exchange densities are comparable to the brane energy density. By 
adopting a phenomenological fluid ansatz and generalizations of it, a set of exact solutions 
of the Friedmann equation that exhibit accelerated expansion were derived. The braneworld 
models with bulk-brane energy exchange allow for crossing of the w = —1 phantom divide 
line without introducing phantom energy with quantum instabilities 26[. The class of the 
braneworld models with bulk-brane energy exchange have the ability to provide crossing of 
the phantom divide, and the observational cosmological data hint towards natural values for 
the model parameters. 

It is the purpose of the present paper to construct a model of braneworld reheating 
by adopting as a key feature the energy exchange between the bulk and the brane during 
the reheating stage. The era of focus will be from the end of inflation until the present 



time, in contrast with 



and 



25 



261 . In order to explain the present day composition 



of our universe, we adopt, as an essential characteristic of our model, the hypothesis that 
the existence of the dark matter is the result of the real matter exchange between the bulk 
and the brane. However, the other matter fields bare no such freedom, according to the 
confinement conditions introduced in j^. 
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The main physical features of our model can be presented as follows: the dark matter 
from the bulk and flows to the brane after the inflationary phase, while the scalar fields and 
the normal matter are always confined on the brane. Therefore the bulk leaves an imprint 
on the brane by injecting dark matter. The effect of dark matter flowing into the brane 
is studied in the context of the cosmological composition, with the reheating temperature 
as an important parameter. The inflaton field decays into the normal matter field [ll|, 2?! . 
attaining thermal equilibrium and acquiring plasma masses through the mechanism proposed 
in 15|. Therefore, as reheating is concerned, the presence of the scalar field is essential, and 



it decays into visible matter only. Furthermore, the specific forms of the inflaton potential 
are not essential in our approach, since during reheating, the kinetic term cfp plays a crucial 
role. The main focus of the present paper is to study the dynamics of the universe by taking 
into account the extra-dimensional effects on the reheating process. The ratio u of the 
energy densities of the dark matter and of the normal matter is an important parameter of 
our model, and its time evolution is studied by using both analytical and numerical methods. 
As a test of the model, several important observational parameters like the maximum density 
of the matter in the universe and the reheating temperature Treh are also estimated. Hence, 
we adopt as our main observational constraints the constancy of the ratio u after the photon 
decoupling, the absence of t 



from gravitino product. 



28 



[he inflaton in the present day universe, and the bound of T^eh 

The present paper is organized as follows. The basic equations of our model are formu- 
lated in Section II. A dimensionless form of the field equations is presented in Section III, 
where the physically acceptable range of the values of the physical parameters is also dis- 
cussed. The extra-dimensional effects dominated phase of the reheating, during which the 
quadratic terms in the energy density dominate the dynamics and evolution of the universe, 
is discussed in Section IV. The reheating process is considered in its full generality in Section 
V, and the role of the linear terms in the energy density is investigated. Finally, in Sec. VI 
we summarize and conclude our results. 
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II. GEOMETRY AND FIELD EQUATIONS IN THE BULK-BRANE ENERGY 
TRANSFER MODEL 



We start by considering a five dimensional {5D) spacetime (the bulk), with a single four- 
dimensional (4D) brane, on which usual (baryonic) matter and physical fields are confined. 
The AD brane world (^^^M,g^^) is located at a hypersurface (i? {X"^) = O) in the 5D bulk 
spacetime C^^^M, qab), of which coordinates are described by X^, A = 0,1, 4. The induced 
AD coordinates on the brane are x'^, fi =J}, 1, 2, 3. 

The action of the system is given by 



S — Sfiulk + Sfyj-anei (l) 



where 



and 



Sbuik = / 

J(5)Af 

Sbrane = / \/-^ 



d^X, (2) 



d^X, (3) 



where fc| = SttGs is the hD gravitational constant, '^^^R and '•^^Lm are the hD scalar curvature 
and the matter Lagrangian in the bulk, Lhrane {dafs, <P) is the 4d Lagrangian, which is given 
by a generic functional of the brane metric Qa/B and of the matter fields 0, is the trace of 
the extrinsic curvature on either side of the brane, and A5 and A (the constant brane tension) 
are the negative vacuum energy densities in the bulk and on the brane, respectively. The 
brane tension is the characteristic energy scale of the brane, below which the quadratic terms 
in the field equations, derived from the action Eq. ([T]), can be neglected. Thus, in the limit 
A ^ 00 the field equations reduce to the usual form of the standard general relativity. 
The Einstein field equations in the bulk are given by 

^'^Gjj = kf'^Tjj, ('^Tjj = -A.^'^gjj + 6(3) [-X^'^gjj + Tjj] , (4) 

where 

- -2^^^^^ + '='9,y^'L„, (5) 

is the energy-momentum tensor of bulk matter fields, while T^i, is the energy-momentum 
tensor localized on the brane and which is defined by 

Till/ — 2 - — -|- Quu Lijrajie- (6) 
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The delta function S (B) denotes the locahzation of brane contribution. In the 5D space- 
time a brane is a fixed point of the Z2 symmetry. In the following we assume 



(7) 



The induced AD metric is gu = ^^^gij — rijuj, where nj is the space-like unit vector field 
normal to the brane hypersurface ^'^^M. The basic equations on the brane are obtained by 



projections onto the brane world, and are given by 



3l| 



G 



A(5) 



1 



1 



where 



S, 



^^lu — '^ABCDn n g^g^, 



and 



F -(5) 



Tab9X + (^'^Tabu 



1 (5) 

^^^-4 T]g,^, 



(8) 

(9) 
(10) 

(11) 



respectively. 

Apart from the terms quadratic in the brane energy-momentum tensor, in the field equa- 
tions on the brane there are two supplementary terms, corresponding to the projection of 
the 5D Weyl tensor e^u and of the projected tensor F^,^, which contains the bulk matter 
contribution. Both terms induce bulk effects on the brane. 

For cosmological applications we adopt a metric of the form 



(12) 



where 7^,^ is the maximally symmetric 3-space, and we denoted the fifth coordinate by 
= S,- a and n are the scale and shift factors, respectively. The Hubble parameter H on 
the brane, describing the cosmological dynamics of the universe, is defined as H = d/ a. For 



further discussions on the gravitational 



cosmological applications see 



32 



33 



ield equations in the brane world models and their 



34 



35 



36 



37| 



On the brane the matter, consisting of normal matter, dark matter, a scalar field etc., is 
in a perfect fluid form, and is described in terms of a total energy density ptotai = XliLiP*' 



where pi 



l,...,n are the energy densities of the individual components, and a total 
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thermodynamic pressure ptotai = ^^=iPi, where pi, i = l,...,n are the pressures of the 
matter fluid components. The corresponding energy-momentum tensor on the brane is 



Tl^ = (5(0diag {-ptotai, Ptotai, Ptotai, Ptotai) ■ (13) 

The delta function implies the confinement of the matter fields on the brane at ^ = 0. 
We assume that the matter component in the bulk is in the form of dark matter, with 
energy density p^M and pressure Pdm, respectively, with energy- momentum tensor 

Ti = {pDAi + Pdm) U^Ub + PdmS^. (14) 

Then the energy momentum tensor of dark matter in the bulk has the non-zero 
components : 

^5° ~ {Pdm + Pdm)U^ , Tl = {pDM + Pdm)U^U^ + Pdm, (15) 



where the five- velocity represents the matter flow from the bulk to the brane [24|. Hence 
in our model the dark matter can exist simultaneously in the bulk and on the brane. 

The term 2T^ is the discontinuity of the (0,5) component of energy momentum tensor at 
^ = 0. In the static bulk with respect to the expanding brane, f/5 is given by oc —IH, 



where / = a/— 6M5 /A5 is the bulk curvature radius, and M5 and A5 are the Planck mass 



38, 



39| . We assume 



and the cosmological constant in the 5 dimensional bulk, respectively 
that the dark matter obeys an equation of state of the form 

Pdm = {w- l)pDM, (16) 

where 1 < w < 2 is a constant. 

n 3^ 

To obtain T5 we follow the approach developed in [24]. Since Pdm = per {c^o/d) where 
per is the present day critical density and qq is the value of the scale factor at the moment 
when the energy transfer from the bulk to the brane begins, we obtain 

T,=--p..(-j H. (17) 

The dimensionless parameter abb describes the flow of dark matter from bulk to brane. 



Finally, the evolution equation of the dark matter on the brane is 



24| 



PDM + SH{pDM + Pdm) = -2T° (18) 



IG 



In the widely accepted inflationary scenario it is assumed that during an initial period 
the universe is dominated by a large, approximately constant potential term V{(j)) of a scalar 
field 0, known as the infiaton field [9|, |lO| . The energy- momentum tensor of the scalar field 
can be written in the perfect fiuid form, with energy density and pressure given by 

p^ = Y + m, (19) 

and 

P0 = y-n0), (20) 

respectively. During a second period of evolution, the potential minimum is approached, 
V{(f)) tends to zero, and the scalar field starts to fiuctuate violently around the minimum 
value. The coupling to the field component becomes effective and the scalar field e nerg y is 
converted into matter energy, thus giving rise to a large increase of the entropy 

Infiationary scenarios have always to face the problem of the transition from a de Sitter 
stage, during which the evolution of the universe is dominated by the scalar field, to a sub- 
sequent radiation - or matter - dominated Friedmann- Robert son- Walker type cosmological 
model. One of the possible approaches to this problem is phenomenological 10, lUl, l27|] , and 
consists in the introduction of a suitable loss term in the scalar field dynamics, which also 
appears as a source term for the energy density of the matter fiuid. It is this source term 
which is responsible for the reheating process that follows the adiabatic supercooling process 
during the de Sitter phase. To describe mathematically the process is convenient to intro- 
duce a "friction" term F, describing the decay rate of the scalar field, and a corresponding 
term for the matter fiuid. 

Therefore, the basic equations describing the reheating phase after braneworld infiation, 
with matter exchange between the brane and the bulk, are given by 

ff^=f^V = 5^(l + ^)-4 + A. + X. (21) 



pm + ^H{pm+Pm) ipm,P^) P<j,, (22) 

p^ + 3H{p^ + p^) = -T {p^,p^) p^, (23) 

abb fdo^.^^ 

Pdm + 3H{pdm + Pdm) = -^Pcr y—J H, (24) 

X + 4^ (x + f T|) = I (p,„,., + A) T°, (25) 



Ptotal = Pm + P4> + PDM, (26) 

where A4 is the cosmological constant on the brane, x is the so-called dark radiation term, 
and we have assumed the existence of a single matter component on the brane, with energy 
density pm and pressure Pm, respectively. 



III. DIMENSIONLESS FORM OF THE FIELD EQUATIONS 

The system of equations Eqs. fl?T]) - fl2B]) can be simplified by adopting some, physically 
justified, approximations. First of all, the AD cosmological constant A4 can be neglected, 
as its derived value from the 3rd year WMAP data is around 10^^'^GeV*^ js^, which is a 
small fraction of the energy density up to at least the photon decoupling period. Also the 
curvature term can be taken as zero. A; = 0, as the universe is observed to be nearly flat. 
Moreover, under the conditions considered in this paper, the dark radiation term 

X oc ^ (27) 



2^. 



is negligibly small as compared with the other components 

An important relation to characterize the matter is the equation of state (EOS) pm = 
(7 — 1) Pm- The product particles have different equations of state in various time scales. The 
index 7 is 7 = 4/3 for radiation before the radiation- matter equality time at teg = lO^years 
~ 10^^ sec, and 7 1 for non-relativistic matter after teq- The constant w determining 
the EOS of dark matter will be fixed in various scenarios. It takes the value w = A/3 for 
hot dark matter (HDM), and w ^ 1 for cold dark matter (CDM), just like in the case of 
ordinary matter. 

Next we consider the scalar field properties. In the post-inflationary stage, the infiaton 



field executes coherent oscillations about the minimum of the potential 4l(|. Therefore the 



kinetic term dominates the potential term in the reheating era, thus leading to 

P0 = y + m^ Y^P^- (28) 



The explicit expression of the decay width of the scalar field can be represented as [Ji 

T = a^MJl-(^] (29) 
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where a^f, and are the couphng constant and the mass of inflaton respectively. Assum- 
ing that the thermodynamic equihbrium is estabhshed, T is the temperature of the decay 
product, and it can be related to the matter density p^- Generally, 

p„ = aT^/(^-i), (30) 

where a is a constant. If the matter is in the form of radiation, pm = 7r^T^/15, and 
cr = 7r^/15. In the matter-dominated universe, the relation between p^ and T can also 
be written down explicitly. The scalar field is negligible small in the matter-dominated 
phase. In the non-relativistic matter domination era T < leV. This is far smaller than 

n 

the minimum bound obtained for [15]. Hence, in the non-relativistic phases of matter 
evolution the decay rate is simply F = a^M^. 
Thus, by defining the parameter 

and with the use of Eqs. ( 1271) - ( l30i) . the system of equations describing the reheating process 
on the brane in the presence of bulk-brane energy exchange is given by 

Ida ( StiG 2 2 \^ /nr.\ 

-^ = H=[^Ptota^ + ^Pto.al (32) 



Pdm = H 




I I/O \2(7-l)/7 

-3^Hp^ + 2a^M^p^Jl-j^[^^) , (33) 



'Pm\2(7-l)/7 



-QHp^-2a^M^p^Jl--^{^) , (34) 



abb faoY^ 



(35) 



In order to simplify the study of the evolution equations during the reheating period, we 
rescale the time variable and the densities, by introducing a set of dimensionless variables 
tdm, r^, Vcr), defined as 

t = Bt, Pm = Arm, Pdm = Atdm, P4, = Ar^, p^r = Avcr, (36) 

where 

R = 



and 

A = (38) 



respectively. By denoting 

the system of equations describing the reheating phase with bulk-brane energy exchange is 
given by 

= ^(^) = \/^m + rDM + r^ + 7]'^ (r„ + Tdm + r^f, (40) 



^ = -37/i(r)r^ + r ^sJl^^n^F^', (41) 



^ = -6/1 (r) - ^1 - o^r^^^"^)/^. (42) 
ar 



/i(r) 



(43) 



The system of equations Eqs. fH0l) - fj42l) must be considered together with the initial 
conditions 

a(0) = ao, rDM(O) = 0, r™(0) = 0, r^(0) = r<^o, (44) 

where ao and r(^o are the initial values of the scale factor and of the scalar field energy density 
at the beginning of the reheating phase. 

In order to obtain a consistent physical interpretation of our results, we need to know 
the numerical values of the model parameters, which are a<^, and A, respectively, aif, 
and are not known separately, but the whole quantity a^M^ can be estimated. From 
Eq.( l29i) . it follows that the order of magnitude of aipM^ is about the order of magnitude of 
the decay width, which is the reciprocal of the characteristic timescale of reheating. The 
inflationary era ends, and reheating can start at the earliest at around t = 10^^^ s, while 



the Hot Big Bang commences at around t = 10 s 



m 



4l| . The reheating process should 
complete before the Hot Big Bang to restore the BBN. Therefore 10^^ sec.~^ < T < 10^^ 
sec.~^ implies that 

1 keV < a^M^ < 10^ GeV. (45) 

As for the brane tension A, it is constrained by the validity of Newtonian gravity in 4 
dimensions on length scales smaller than 0.1mm 4^. The minimum bound is given by 



A > (100 GeVy. (46) 
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With the use of the second transformation in Eqs. (136|) we can express the dimensionless 
matter density for an ultra-relativistic gas (7 = 4/3) and with energy- density pm = 
7r^T^/15 as a function of the temperature T as = An^GT'^/dOa'^M^. From the definition 
of the decay width we obtain the constraint T < on the temperature. Therefore, for a 
universe in which the normal matter is in the form of radiation only, this condition can be 
reformulated in a form of a constraint on the dimensionless matter density r^, given by 

< -mt- '-''^ 

On the other hand from Eqs. (HTl) and (jl2l) we obtain the condition curm < 1, which is 



consistent with Eq. (1471) . 



IV. EXTRA-DIMENSIONAL EFFECTS DOMINATED BRANE REHEATING 

At the end of the inflation, there is no ordinary matter present in the universe. Even in 
the first early stages of the reheating, the scalar field dominates the matter energy densities, 
P0 ^ Pm,PDM- Due to the high energy of the scalar field, the extra-dimensional effects 
due to the presence of the fifth dimension are important at this stage, and the quadratic 
terms in the evolution equations dominate the cosmological dynamics. Therefore one can 
neglect the linear term in Eq. ( HOl) with respect to the quadratic term. Hence, the extra- 
dimensional effects are dominant during the reheating era. Moreover, we assume that the 
condition T << is satisfied during the entire reheating process. This condition implies 
ujrm^~^^^"' << 1. Then the equations describing the reheating process with bulk-brane 
energy transfer are given by 

-3- = ^ (^m + roM + r^) , (48) 
a ar 

droM , „ Ida (ao\^^lda . . 

—, — ^^w-—rDM = e\ — ] - — , (51) 
dr a ar \ a / aar 

Eq. ( |50l) can be immediately integrated, and with the use of the initial conditions we 
obtain the evolution of the energy density of the scalar field during reheating as 

r,(r)=r,o(^)%-^ (52) 
13 



Eq. (15T!) gives the density of the dark matter on the brane in the form 
The evolution of the matter energy density may be represented as 



(53) 



(54) 



A. The stiff matter case 



The value 7 = 2 corresponds to the case of the so-called stiff matter fluid, with equation 
of state given by pm = pm- In this case the energy density of the matter is given by 

- r,o(l-e--)(^)',7 = 2. (55) 



The sum of the energy densities of the matter and scalar field satisfies the equation 



For 7 = 2 we obtain 



d , , 



(37^n 



„ \da 
6r^)- — . 
a dr 



(56) 
(57) 



The scale factor of the universe can be obtained from the general equation 
1 da 
a dr 



^37-6 



e '^dr — P (flo) 



+ 



ao 



a 



(58) 



In the case in which both the dark matter and the usual matter are in the ultra-relativistic 
state, 7 = M7 = 2, and the equation describing the scale factor dynamics is given by 



■ da 
dr 



r^o + e\n 



ao 



,7 



w 



(59) 



In this case the time evolution of the universe during the extra-dimensions dominated 
phase of the reheating era can be expressed as 



exp 



6r, 



00 



Ei ^ + 61n - , (60) 
er] ' \ e J ye V'^o/J 

where Ei (2;) = — j^^e~*dt/t is the exponential integral function, and where the principal 
value of the integral must be taken. In a power series representation the dimensionless time 
parameter r can be obtained in terms of the scale factor a as 



1 fa 
ei] \ao 



1 



Qr^o/e + 61n (a/ag) 



+ 



1 



Qr^o/e + 61n(a/ao) 



+ ... 



(61) 
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B. The radiation dominated phase 

In order to obtain a clear physical picture of the reheating process during the extra- 
dimensional effects dominated era one must study the evolution equations Eqs. (jUD - ([51]) 
numerically. In order to do this we assume that the normal matter is always generated from 
the inflaton field in the form of ultra-relativistic particles, leading to a very hot and high 
density initial state of the universe, an assumption which is consistent with the standard 
hot Big-Bang model. Therefore we take 7 = 4/3. As for the dark matter, it can be injected 
from the bulk both as a relativistic fluid, with t/; = 4/3, or as a non-relativistic, pressureless 
fluid, with w = 1. The variation of the radiation energy density is represented, for both 
equations of state of the dark matter, for fixed values of the parameters e and u and for 
different values of the parameter rj in Fig. 1. 




FIG. 1: Variation of the energy density of the radiation (7 = 4/3) during the extra-dimensional 
effects dominated reheating period, with ultra-relativistic dark matter injection from the bulk 
{w = 4/3) (figure on the left) and non-relativistic pressureless dark matter injection (w = 1) (figure 
on the right), respectively, for to = 0.5, e = 6 and different values of r/: rj = 0.5 (solid curve), rj = 1 
(dotted curve), = 1.5 (short dashed curve) and rj = 2 (long dashed curve). The initial values of 
the scalar field density and of the scale factor are r^Q = 1 and 00 = 10~^, respectively. 

The variation of the dark matter energy density is represented in Fig. 2. 

The variations of the scalar field energy densities and of the scale factors in these two 
situations are represented in Figs. 3 and 4, respectively. 

As one can see from the Figs. 1 and 2, during the reheating period the energy of the 
scalar field is transferred to the normal matter (radiation). The energy density of the 
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FIG. 2: Variation of the energy density of the ultra-relativistic dark matter {w = 4/3) (figure on the 
left) and of the pressureless non-relativistic dark matter {w = 1) (figure on the right), respectively, 
during the extra-dimensional effects dominated reheating period, for lu = 0.5, e = 6 and different 
values of 77: r/ = 0.5 (solid curve), ry = 1 (dotted curve), rj = 1.5 (short dashed curve) and rj = 2 
(long dashed curve) . The initial values of the scalar field density and of the scale factor are r^o = 1 
and ao = 10~^, respectively. The normal matter is in the form of radiation with 7 = 4/3. 




0.2 0.4 0.6 0.8 1 1.2 1.4 0.2 0.4 0.6 0.8 1 1.2 1.4 

FIG. 3: Time variation of the energy density of the inflaton field during the extra-dimensional effects 
dominated reheating period, with ultra-relativistic dark matter injection from the bulk (w = 4/3) 
(figure on the left) and non-relativistic pressureless dark matter injection {w = 1) (figure on the 
right), respectively, for co = 0.5, £ = 6 and different values of ry: t] = 0.5 (solid curve), 77 = 1 (dotted 
curve), rj = 1.5 (short dashed curve) and r/ = 2 (long dashed curve). The initial values of the scalar 
field density and of the scale factor are r^o = 1 and ao = 10~^, respectively. The normal matter is 
in the form of radiation with 7 = 4/3. 

matter increases from zero to a maximum value, and, after reaching the maximum, it starts 
to decrease, due to the cosmological expansion. During the same period there is an injection 
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FIG. 4: Time variation of the scale factor of the universe during the extra-dimensional effects 
dominated reheating period, with ultra-relativistic dark matter injection from the bulk (to = 4/3) 
(figure on the left) and non-relativistic pressureless dark matter injection {w = 1) (figure on the 
right), respectively, for uj = 0.5, e = 6 and different values of ry: rj = 0.5 (solid curve), rj = 1 (dotted 
curve), r/ = 1.5 (short dashed curve) and rj = 2 (long dashed curve). The initial values of the scalar 
field density and of the scale factor are r^Q = 1 and oq = 10"'^, respectively. The normal matter is 
in the form of radiation with 7 = 4/3. 

of dark matter from the bulk, the energy density of the dark matter increasing from an 
initial zero value to a maximum value, which depends on the equation of state of the dark 
matter and on the parameter rj, describing the magnitude of the extra-dimensional effects. 
During the reheating process the energy density of the scalar field decreases to zero. 

V. GENERAL EVOLUTION OF THE BRANE UNIVERSE DURING THE RE- 
HEATING PERIOD 

After reaching its maximum value during the reheating period, due to the cosmological 
expansion, the energy density of the matter (in all its forms) decreases in time and for 
time intervals t > tcr the condition ptotai/^ << 1 is satisfied. In this case the quadratic 
contributions to the field equations can be neglected, and the universe is described by the 
standard general relativistic cosmological models. In the dimensionless formulation of the 
field equations the general relativistic limit corresponds to r] 00. In the general case, 
when all terms are considered in the field equations, the evolutions of the matter and of the 
scalar field is given generally by Eqs. fl32]) - (jM]), respectively. However, the linear term in 
energy density modifies the expansion rate of the universe, which is described by Eq. (I40l) . 
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From the mathematical structure of the field equations one can derive some general prop- 
erties of the model, and of the cosmological evolution during the reheating. The normal 
matter and the dark matter reach a maximum value at a moment r = Tmax, the maxi- 
mum values of the densities being given by the solutions of the equations drm/dr = and 
f^DAi/d'^ = 0, respectively. In the case of the dark matter the maximum value is given by 

e / ao ^ 



r 



(max) 
DM 



3w \ a 



(62) 



The maximum value of the dark matter density is determined by the injection rate e, the 
cosmological dynamics, determining value of the scale factor at Xmax, and by the equation 
of state of the dark matter. By assuming that c^r^^"'^/^ « 1, the maximum value of the 
normal matter density can be obtained from the algebraic equation 



(max) 



87/1 (r) 

In the case of the extra-dimensional effects dominated evolution, 
V {.fm + ruM + r^) , and Eq. ([63]) gives 
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which can be approximated as 
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(64) 



(65) 



The maximum value of the matter energy density is inverse proportional to the parameter 
rj = a/3/ AnGXa^M^p. In the case of a fixed value of a^M^, a large value of rj, corresponding, 
for example, to a small brane tension A, would determine a decrease of the maximum matter 
density, while small values of 1] (large brane tensions) would increase the amount of matter 
produced during the reheating phase. In terms of the scale factor the maximum matter 
density can be expressed as 



r 



(max) 



1 



ao 



3«)— 6 



0,-0 . ^ - . - (66) 

3777 I owr^ ^ 1^ \ I 

Generally, the ratio u of the dark matter and normal matter energy density can be 
represented, as a function of the scale factor, as 

rDM jr) ^ ^ 37-6 /aoy"'-37 In (a/ao) 



u 



In the small time limit the energy density of the baryonic and of the dark matter is 
negligibly small with respect to the energy density of the scalar field, which in the very early 
stages of the reheating is given by Eq. fl52p . Therefore in the limit of r << 1, the evolution 
of the scale factor is obtained as a ~ Oq (6?7r^o)^^^ t^^^- Hence the small time behavior of u 
is described by the equation 

^'00 

Since in the limit of large time the function a^'^~^e~^ tends to zero, it follows that the 
value of the integral a^'^'^e^'^dr can be written as a function of the initial value of the 
scale factor ao, so that /J" a^'^^^e^^dr ^ /5 (ao). 

Therefore in the large time limit the ratio of the densities of the dark matter and of the 
baryonic matter can be written as 

eaf?~'^ /ao\3"'"37 / a \ , , 

u = ^-^l— ) In — ,r>>l. (69) 



r<poP (ao) ^ a / \aQ 
Precision cosmological observations have shown that u has a constant numerical value, 
at a high degree of accuracy. Since in the present day universe the baryonic matter is 
non-relativistic, with 7 = 1, the constancy of u requires for the dark matter to also be in 
a non-relativistic (cold) phase, with w = 1. The magnitude of the dark matter /baryonic 
matter ratio depends on the initial conditions of the universe during the reheating period, 
namely, the dark matter injection rate, the energy density of the scalar field and the value 
of the scale factor. On the other hand, the ratio of the densities of the dark matter and of 
the normal matter is not an absolute constant during the cosmological evolution, but it is a 
very slowly, logarithmically increasing function of time. 

The results of the numerical integration of the Eqs. ( HOl) - ( H3l) are represented in Figs. 
5-8, respectively. Fig. 5 presents the variation of the normal matter energy density during 
reheating, with the linear terms in the density also included, for the two equations of state 
of the dark matter. The variation of the dark matter energy density is represented in Fig. 
6. 

The variations of the scalar field energy densities and of the scale factors for the two 
equations of state are represented in Figs. 7 and 8, respectively. 

As one can see from the figures, the inclusion of the linear matter terms into the evolution 
equations leads to an overall decrease of the numerical values of the cosmological parame- 
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FIG. 5: Variation of the energy density of the photon gas (7 = 4/3) during the reheating period, 
with ultra-relativistic dark matter injection from the bulk {w = 4/3) (figure on the left) and 
non-relativistic pressureless dark matter injection [w = 1) (figure on the right), respectively, for 
00 = 0.5, e = 6 and different values of 77: r] = 0.5 (solid curve), 77 = 1 (dotted curve), 77 = 1.5 (short 
dashed curve) and 77 = 2 (long dashed curve). The initial values of the scalar field density and of 
the scale factor are r^o = 1 and oq = 10~^, respectively. 
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FIG. 6: Variation of the energy density of the ultra-relativistic dark matter {w = 4/3) (figure on 
the left) and of the pressureless non-relativistic dark matter {w = 1) (figure on the right) during 
the reheating period, for lu = 0.5, e = 6 and diflFerent values of 77: 77 = 0.5 (solid curve), 77 = 1 
(dotted curve), rj = 1.5 (short dashed curve) and 77 = 2 (long dashed curve). The initial values of 
the scalar field density and of the scale factor are r^o = 1 and oq = 10~^, respectively. The normal 
matter is in the form of radiation with 7 = 4/3. 

ters during the reheating period. If in the case of the extra-dimensional effects dominated 
dynamics the maximum values of the matter energy density are of the orders of ~ 0.12, 
in the general case and for the same value of the physical parameters this value is ~ 0.06. 
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FIG. 7: Time variation of the energy density of the inflaton field during the reheating period, 
with ultra-relativistic dark matter injection from the bulk {w = 4/3) (figure on the left) and 
non-relativistic pressureless dark matter injection (w = 1) (figure on the right), respectively, for 
uj = 0.5, e = 6 and different values of 77: rj = 0.5 (solid curve), r] = 1 (dotted curve), r] = 1.5 (short 
dashed curve) and rj = 2 (long dashed curve). The initial values of the scalar field density and 
of the scale factor are = 1 and oq = 10~^, respectively. The normal matter is in the form of 
radiation with 7 = 4/3. 
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FIG. 8: Time variation of the scale factor of the universe during the reheating period, with ultra- 
relativistic dark matter injection from the bulk {w = 4/3) (figure on the left) and non-relativistic 
pressureless dark matter injection {w = 1) (figure on the right), respectively, for uj = 0.5, £ = 6 and 
different values of 77: r) = 0.5 (solid curve), r) = 1 (dotted curve), 77 = 1.5 (short dashed curve) and 
r) = 2 (long dashed curve). The initial values of the scalar field density and of the scale factor are 
r^o = 1 and ao = 10~^, respectively. The normal matter is in the form of radiation with 7 = 4/3. 



However, the maximum value of the energy density of the dark matter is not significantly 
influenced by the change in the cosmological dynamics, since it is essentially determined by 
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the injection process from the bulk. On the other hand the decay of the inflaton field is 
more rapid when the full dynamics on the brane is taken into account. 

It is easy to show that the ratio of the dark matter and normal matter energy density u 
satisfies the following evolution equation 

"(3^ _ 3^) hir) - ^Vl-^^^^'^'l - + 4-) " (70) 



du 
dr 



The maximum value of m, Umax — u (rmax) is given by the solution of the equation 
du/dr]^^^ — 0, giving for the maximum value of u the expression 



- (37 - M h{r)rr. 



(71) 



— Tmax 

In the limit of large times, when r<^ 0, and assuming that the dark matter obeys the 
same equation of state as the normal matter, w — ^, u satisfies the equation 

^^sf^V'"^. (72) 
dr \ a / 

For non-relativistic dark matter and normal matter we obtain the simple relation 

du , , . 

— =£/i(T). (73) 

Thus the variation of the ratio of the dark matter and of the normal matter is proportional 
to the Hubble parameter. 

To study the time variation of the ratio of the dark and baryonic matter u we consider 
that the transition between the relativistic and non-relativistic state can be modelled as a 
smooth transition, with the equations of state of the form pdm — w{t)pDM and Pm — li't)Pm, 
respectively. In the dimensionless time variable r we assume that the functions 7(r) and 
w{t) can be given as 

or + To J T + To-y 

and 

3 r + To^ r + row 

respectively, where and Tow are the critical transitions times for the dark matter and the 
baryonic matter. When r << To-y,ro^, 7 ~ 4/3, while for r >> To^,Tou,, 7 ~ 1. The time 
variation of u is represented, for two different numerical values of the critical transitions 
times of the dark matter and baryonic matter, respectively, and for different values of the 
physical parameters, in Fig. 9. 
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FIG. 9: Time variation of the ratio u of the dark matter and normal matter energy densities during 
the reheating period, by assuming a smooth transition from the relativistic to the non-relativistic 
state, so that 7(r) = (4/3)ro^/(T + Toy) + r/(r + tq^) and w{t) = (4/3)tou,/(t + tow)t/{t + tqw), 
respectively. The values of the parameters are co = 0.5, rj = 2 (solid curve), rj = 2.5 (dotted curve), 
77 = 3 (short dashed curve) and = 3.5 (long dashed curve). We have assumed that the transitions 
from a non-relativistic to a relativistic phase takes place at the critical time tq^ = tqw = 10^. 
The initial values of the scalar field density and of the scale factor are r<^o = 1 and ao = 10~^, 
respectively, while the dark matter injection rate is e = 0.01. 

VI. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have considered the reheating of the universe in the brane world 
models with bulk-brane energy exchange. Our model is essentially based on the hypothesis 
that the present day dark matter, representing around 25% of the total matter content of 
the universe, was transferred onto the brane during the post-inflationary reheating period, 
while the normal matter is the result of the decay of the infiaton field into a ultra-relativistic 
photon gas. Hence, after the inflationary epoch nor dark matter neither normal matter were 
present on the brane, and they were created during the reheating. 

To study the reheating process we have introduced a set of dimcnsionless quantities, which 
describe the evolution of the scaled densities in terms of a dimcnsionless time parameter r. 
By using the results of the numerical simulations for the rescaled variables one can obtain 
some constraints on the physical parameters of the model. 

The reheating time trehi the time after which the universe is reheated to a temperature 
allowing the standard hot big bang model to take place is related to the dimcnsionless time 
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Treh by Eqs. (l36i) and (1371) . respectively. These equations allows us to constrain the mass 
and the coupling constant of the infiaton field as 

a^M^ = ^ = 0.328Trehtreh{s) X lO'^^ GcV. (76) 

Since according the hot standard big-bang scenario reheating must end at around treh = 
10^^^ s, it follows that a^M^ = 0.328Treh keV. In our model and for the typical values 
of the physical parameters we have adopted (bulk-brane injection rate and brane tension), 
according to the numerical analysis reheating may definitely take place for values of Treh 
in the range 0.1 < Treh < 1 (see Figs. 1-2 and 5-6), thus giving 0.0328 keV < a^M^ < 
0.328 keV. Values of a^M^f) of the order of 1 keV or higher would require values of Treh of 
the order of Treh > 3. 

The numerical analysis of the field equations, as well as their qualitative study, shows the 
existence of an absolute maximum value of the matter density p^mattlr during the cosmological 
evolution of the universe. The maximum value of the total matter density is given by 



(max) 
Pmatter 



(max) I (max) 
' m "T ' DM 



SttG 



7.13 [r^'^^) +r™J (a^M^/GeV)' x 10^^ GeV^ 

(77) 

By assuming that r^^^^ + r^^^^^ ~ 0.5 — 0.7 and Q;(/,M(^ ^ 1 keV, we obtain 

3.565 X 10^5 GeV^ < p^^r < 4-991 x 10^^ GeV^ (78) 

Eq. (ITHj) allows the determination of the reheating temperature Treh of the universe, once 
the equation of state of the matter is known. Assuming that both normal matter and dark 
matter are in a ultra-relativistic state, we obtain 

Treh ^ 3 X 10*^ GeV. (79) 

Several estimations of the numerical value of the reheating temperature have been ob- 
tained by using some observational constraints. A strong constraint on the cosmological 
models comes from the study of the CMB anisotropy, and requires that the results do not 



exceed the observed level of normalized anisotropy on large scales |ld . Il5[ |. By also re- 
quiring that the radiative corrections mediated by the couplings do not spoil the fiatness 
of the potential limits, it follows that the reheating temperature must be below the grand 
unification theory scale, Treh < 10^^ GeV, which physically means that the grand unified 
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theory symmetries are not restored and therefore there is no second phase of production 
of monopoles 1^, [l^. By assuming that the energy density of the universe is in the form 
of relativistic matter with density ptotai = fi'*7r^T^/30, where g^, is the effective number of 
massless degrees of freedom ((7* = 100 — 1000), the reheating temperature can be obtained 
as T,eh ^ Q.2{im/ g^f'^^/Vtotmpi ^ 0.2 [im/ g,f'^ a/^ ^ M^mpu where is the total 
decay width and mpi is the Planck mass. From the study of the CMB anisotropics it follows 
that ^ IQ-^mpi [lo|, ij- 

In the context of supersymmetry models, a stronger constraint may be obtained from the 
production rate of gravitinos, the supersymmetric partners of the gravitons. The gravitinos, 
having a long lifetime, of the order of 10^ s, may cause serious problems in cosmology. For 
a wide range of the gravitino masses in order to obtain consistency with the nucleos ynt hesis 
constraints, the reheating temperature must be below 10^ GeV, Treh < 10^ GeV [lO|, l30 |. 
A stronger constraint on the reheating temperature for the gravitino decay products not to 
destroy the light elements synthesized during the Big-Bang nucleosynthesis was obtained in 



45|, and requires that Treh must satisfy the condition Treh < 10^ — 10^ GeV. Primordial 



nucleosynthesis requires that the universe is close to thermal equilibrium at a temperature 



of around 1 — 4 MeV ll|] . Hence, the reheating temperature predicted by our model satisfies 
all the CMB and nucleosynthesis constraints. 

By adopting the CMB constraint implying that ^ 10~^mpi ^ 1.22 x 10^^ GeV, we 
can estimate the infiaton field coupling constant as a^M^ = 0.328Trehtreh{s) x lO^^'^/M^ = 
2.69 X Trehtreh{s) X 10""^*^. For Ueh ~ lO""^*^ s wc obtaiu ^ 2.69 x Treh X 10~^^, or 
2.69 X 10^2^ <a^< 2.69 x lO^^e^ 

The numerical values of the cosmological parameters sensitively depend on the parameter 
7], a dimensionless combination of the mass and coupling constant of the infiaton field and 
of the brane tension, respectively. From Eq. (!38|) we obtain the brane tension A as 

A = = ^ (a,M,/GeV)^ x 10^^ GeV^ (80) 

Since from the results of the numerical integration it follows that the values of rj are of 
the order of 0.5 < 1] < 2, and by adopting again for a^M^f, a value of 1 keV, we obtain 
0.891 X 10^^ GeV^ < A < 14.26 x 10^^ GeV^. This value is much larger than the minimal 



value A > 10^ GeV previously obtained in [43|, |4^, but it is consistent with the bound 
A < 10^0 GeV^ derived in |46|. 
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In the present model the reheating temperature satisfies with a very good approximation 
the relation 

Treh = (81) 

relating the reheating temperature to the brane tension, and which was first suggested in 
47|. The universe exits the extra-dimensions dominated phase, when ptotai » A, at a 
density ptotai ~ A, which, assuming that the matter and field content is formed from ultra- 
relativistic particles, gives a transition temperature Ttrans ~ A^^^. According to our results, 
Ttrans ~ Treh- Therefore, the reheating phase of the universe essentially took place during 
the extra-dimensions dominated cosmological era, when the bulk effects may have played a 
major role in the dynamics and cosmological defrosting of the universe. On the other hand, 
since the reheating temperature may be estimated from observations, the exact relation 
between T^eh and A gives the possibility of the direct determination of the brane tension 
from cosmological observations. 

The basic feature of the present reheating model is the continuous energy-matter transfer 
between the brane and the bulk. Even that this transfer was maximum during the reheating 
period, essentially contributing in establishing the present day cosmological chemical com- 
position of the universe, a very small rate of change of the ratio of the dark matter and of 
the normal matter still occurred during the later phases of the cosmological evolution. The 
present day rate of change of this ratio can be estimated from Eq. (173|l . and is given by 

^^ex Ho^ 3.2Ae x 10"^^ s'S (82) 

where Hq is the present day value of the Hubble parameter. Since in order to explain the 
cosmological ratio of the matter components e should have numerical values of the order 
of unity, it follows that the time changes of the dark matter and normal matter ratio are 
negligible small in the recent cosmological times. 
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